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We investigate the isospin-violating mixing of the light scalar mesons ao(980) and /o(980) within 
the unitarized chiral approach. Isospin-violating effects are considered to leading order in the quark 
mass differences and electromagnetism. In this approach both mesons are generated through meson- 
meson dynamics. Our results provide a description of the mixing phenomenon within a framework 
consistent with chiral symmetry and unitarity, where these resonances are not predominantly qq 
states. Amongst the possible experimental signals, we discuss observable consequences for the 
reaction J/ 1 !/ — * (f>Tv°n in detail. In particular we demonstrate that the effect of ao — fo mixing is by 
far the most important isospin-breaking effect in the resonance region and can indeed be extracted 
from experiment. 
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I. INTRODUCTION 

Although the light scalar mesons ao(980) and /o(980) 
have been established as resonances long ago, there is 
still a heated debate going on in the literature regarding 
the very nature of these states. Naively one might assign 
them a conventional qq structure, however, at present no 
quark model is capable of describing both states simul- 
taneously as qq states — see, e.g., Ref. [1]. On the other 
hand, as early as 1977 it was stressed that especially in 
the scalar channel the interaction of four-quark systems 
(two quarks, two antiquarks) is attractive [2]. Some au- 
thors have found indications for the existence of compact 
four-quark states [3, 4]. However, the same short-ranged 
interaction can also be the kernel to the scattering of 
pseudoscalars, giving rise to extended four-quark states 
that one might call hadronic molecules or extraordinary 
hadrons [5-8]. Independently, a similar conclusion was 
found in different approaches [9-14]. 

A very different approach to quantify the nature of 
scalar states was presented in Ref. [15], where it was ar- 
gued that the value of the effective coupling constant 
of a resonance to a particular continuum channel is a 
direct measure of its molecular component if the corre- 
sponding threshold is very close to the resonance posi- 
tion. When applied to the case of the /o(980) also this 
model-independent analysis revealed that this scalar is 
(to a very large degree) of molecular nature [15]. This 
picture was further supported by analyses of the reactions 
(j) ->■ 7r°7r°7 and f -> 77 [16, 17]. For the a (980), on the 
other hand, no clear picture emerged from these studies 
(note that also the data is of poorer quality). This might 
either mean that the physical ao(980) has some sizeable 
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FIG. 1: Graphical illustration of the leading contribution to 
the ao — fo mixing matrix element Al defined in Eq. (1). 



admixture of something different from KK or is a virtual 
state. Here more information is urgently called for. 

An observable supposedly of high sensitivity to the 
structure of the scalar mesons was identified at the end 
of the 1970s, when Achasov and coworkers observed that 
the isospin-violating mixing of the isovector 00 and the 
isoscalar fo should be significantly enhanced due to the 
proximity of the kaon thresholds to the poles of both 
mesons. In Ref. [18] it was demonstrated that the lead- 
ing piece of the ao — fo mixing amplitude can be written 
as 

Al = (fo\T\a ) = i g foKR g aoK R Vs (p K ° - Pk+) 

where Pk",+ denotes the modulus of the relative mo- 
menta of the neutral and charged kaon pairs, respectively, 
and the effective coupling constants are defined through 
Frkk = 9rkkP k ' R = a 0' fo- Obviously, this leading 
contribution is just the difference of the unitarity cut con- 
tributions of the diagrams shown in Fig. 1 and is there- 
fore model independent. In addition, the signal is pro- 
portional to the effective coupling constants of the scalar 
mesons that encode the essential structure information, 
as outlined above. As already stressed in Ref. [18], the 
contribution shown in Eq. (1) is unusually enhanced be- 
tween the K + K~ and the K°K° thresholds, a regime of 
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FIG. 2: Graphical illustration of the subleading contribution 
to the ao — fo mixing matrix element Ay defined in Eq. (3). 
The crosses denote isospin-violating vertices. 



FIG. 3: Soft photon exchange contribution to the ao — fo 
mixing amplitude. 



only 8 MeV width. Here it scales as 
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m d 



(2) 



where m u , rrid, and m s denote the current quark masses 
of the up, down, and strange quark, respectively, rh = 
(m u + md)/2, and we neglect electromagnetic effects in 
the kaon masses for this rather symbolic formula. This 
is in contrast to common isospin-violating effects 1 which 
scale as (m<j — m u )/m s , since they have to be analytic in 
the quark masses. It is easy to see that away from the 
kaon thresholds A returns to a value of natural size. In 
the phenomenological calculation of Ref. [19] this effect 
was confirmed. However, also there the different kaon 
masses were the only source of isospin violation. 

A subleading contribution to the mixing amplitude is 
given, in the resonance picture, by isospin-violating cou- 
plings of the resonances to the two-kaon continuum g^ K ^ 
as depicted in Fig. 2, 



Ay = 



(9f KK 9a Q KK + QfoKK 9^ K g)^ PK + 0(p 2 K 



(3) 



Although these effects are regular in the isospin-breaking 
parameter, i.e. of order — m u , they are kinematically 
enhanced due to the unitarity cut oc pk- An assessment 
of the size of such effects obviously relies on an estimate of 
the isospin-violating couplings g^ K j^ that we will provide 
in this paper. It should be stressed, however, that as 
their isospin conserving counterparts, the g^ K ^ are well 
defined, observable quantities. 

In addition there can also be mixing through the ex- 
change of soft photons in meson loops — see Fig. 3 - 
giving rise to the mixing amplitude Ap. The full mixing 
amplitude is then given by 



A = A, 



A. 



A, 



(4) 



If ao and fo had a significant admixture from elementary 
scalars, one should in addition expect a direct ao — fo 
transition to appear. However, this is not included in 
our approach. 



1 By common isospin-breaking effects we refer to those effects that 
occur at the Lagrangian level. 



Our purpose is to improve on the theoretical under- 
standing of the possible ao — fo mixing phenomenon. Es- 
pecially we would like to get a first quantitative under- 
standing of the possible impact of isospin-violating cou- 
plings. To do so we employ the chiral unitary approach, 
in which we will now include isospin-violating effects to 
leading order in chiral perturbation theory (ChPT) in 
both the strong and the electromagnetic sector. This al- 
lows us to address the following issues: 

1. Does the enhancement of isospin violation in the 
effective couplings near the KK thresholds, see 
Eq. (3) , lead to similarly sizeable effects as the kaon 
mass differences? 

2. What is the effect of soft photon exchange on the 
mixing amplitude, see Fig. 3? 

3. What is the resulting mass dependence of the signal 
for the mixing? 

4. In Refs. [20, 21] it was claimed that, in the mass 
range considered, ao — fo mixing is the by far dom- 
inant isospin-violating effect as it emerges from the 
overlap of two narrow resonances with very nearby 
masses. This issue is discussed in Sec. IIIB. As we 
will see, we are now in the position to check this es- 
timate within a dynamical approach for a specific 
reaction. 

For our calculations we use the chiral unitary approach 
as developed in Refs. [14, 22], which provides the am- 
plitudes for the scattering of two pseudoscalars from 
the coupled channel unitarization of the leading-order 
chiral Lagrangian. What is interesting about this ap- 
proach is that it is not only able to describe the data 
on pseudoscalar-pseudoscalar S-wave scattering up to 
1.2 GeV remarkably well, with just one cutoff of natural 
size, but also to dynamically generate the poles associ- 
ated with the lightest scalar mesons without the need to 
introduce them explicitly in the Lagrangian. Hence one 
can avoid a priori assumptions about the nature or even 
the existence of those resonances that come out natu- 
rally as a consequence of chiral symmetry and coupled 
channel unitarity. This reduces the model dependence 
of the approach considerably. In order to establish the 
nature of the generated poles, additional theoretical in- 
formation is necessary, e.g. in Ref. [6] the leading 1/N C 
behavior was used to provide evidence for a non-qq na- 
ture of the scalars. The scattering amplitudes obtained 
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this way have then been used to implement the final- 
state interaction in the next-to-leading-order calculation 
of scalar form factors [23, 24]. 

Another relevant aspect of unitarized ChPT is that it 
can be extended to higher orders, and, indeed, it is also 
possible to use the fully renormalized next-to-leading- 
order unitarized ChPT scattering amplitudes [25-28] to 
match the final-state interactions into form factors [29]. 
However, for the scalar form factors we are interested in, 
it is possible to simplify the approach by matching their 
next-to-leading-order calculations just with the leading- 
order chiral unitary scattering amplitudes. Indeed, it has 
been shown [23, 24] that this approximation already pro- 
vides a very good description of the existing data on 
isospin-conserving processes, and thus provides a well 
founded starting point for our approach. 

In the literature various reactions are discussed that 
should be sensitive to the isospin- violating a — f mixing; 
amongst those are jp — ► pTr°n [30], n~p — > TT nn [31, 32], 
pn -> dTT°n [33-35], dd -> an ?] [36], and J/* -> 
<f)ir°n [37, 38]. In the first three certain differential ob- 
servables are sensitive to a — fo mixing, for the last two 
the cross section is proportional to the square of the 
mixing amplitude, since the corresponding amplitudes 
vanish in the isospin limit. In this paper we will fo- 
cus on the last reaction since the recent measurement 
by the BES collaboration of the isospin-conserving chan- 
nel J — > 4>ir + ir~ [39] shows a very pronounced signal 
of the / . 2 In addition these data were already analyzed 
within the unitarized chiral approach in Ref. [24], go- 
ing back to the formalism developed in Ref. [23], which 
is very convenient for our purposes since it determines 
the isospin-conserving part of our formalism rather ac- 
curately. A variant of this analysis, including isospin- 
breaking sources, forms the basis of our study. 

II. THEORETICAL FRAMEWORK FOR THE 
ISOSPIN-CONSERVING CASE 

A. Isospin states and scalar form factors 

In all reactions listed above, where signals of ao — fo 
mixing are expected, there are three strongly interacting 
particles in the final state. In principle this necessitates 
a full three-body treatment using relativistic Faddeev 
equations. However, since we will focus on a phenomenon 
that occurs within a very narrow kinematic window, we 
will adopt the usual approximation that kinematics can 
be chosen such that the interactions within the particle 
pair of interest can be isolated. This approximation has 



2 In principle one could also study J/ty — > um°r], however, there is 
no clear signal of the fo visible in the corresponding isospin- 
conserving channel J/ty — » lotttt, and therefore also no pro- 
nounced mixing amongst ao and fo should be expected in this 
channel. 



already been demonstrated to provide a good description 
of the data [23, 24] , but should be checked within a Dalitz 
plot analysis of the data, once available. 

Therefore, in what follows we assume that only the in- 
teraction of the two outgoing pseudoscalar particles needs 
to be considered. The full production amplitude is then 
given by the scalar form factor times at most a polyno- 
mial [40]. Since we are interested in a phenomenon that 
occurs within a mass range of a few tens of MeV only, we 
can safely use a polynomial of zcroth order only. 

The considerations above are rather general and should 
hold for all reactions listed above as possible candidates 
to find signals of ao — fo mixing. As a concrete example 
and since it will be used below, we now briefly reiter- 
ate the formalism of Refs. [23, 24] to describe the decay 
of the J/* into a <f) an d two pseudoscalars. We use a 
Lagrangian coupling the two vector particles to scalar 
currents of zero isospin in the form 

C = C^tf"^ [ss + X^nn] , (5) 

where nn — (uu + dd)/V2. The two parameters 
are a priori unknown and have been extracted by fits to 
experimental data in Ref. [23, 24]. Note that these anal- 
yses assume A^ to be real, which amounts to neglecting 
left-hand cuts or crossed-channel final-state interactions 
(the latter can in principle be separated in a careful anal- 
ysis of the Dalitz plot). The matrix element for the full 
decay then involves matrix elements of the scalar cur- 
rents between the vacuum and two pseudoscalars, which 
are described by scalar form factors as follows: 

-V2BT^{s) = (0\nn\nn) I=0 . 
-s/2BT n K {s) = (0\nn\KK) I=o , (6) 
-V2BT n n (s) = (0|fm| W )j =0 , 

and equivalent definitions of the strange scalar form fac- 
tors with the replacements fin — > ss, T™ — > Tf. The 
connection of the two-meson states of definite isospin to 
the basis of physical particles is given in Appendix A. 

B. Form factors and unitarization 

Using ChPT to a certain (e.g. ncxt-to-leading) order 
to calculate the form factors defined in Eq. (6) guaran- 
tees that we obtain a consistent low-energy expansion, 
with the correct chiral loop corrections. However, we 
are interested in the energy region of the ao and fo res- 
onances, which is outside the realm of applicability of 
ChPT amplitudes. The latter are, up to branch cuts 
generated by Goldstone boson dynamics, just polynomi- 
als in energy, and as such cannot generate the poles that 
quantum field theory requires to be associated with res- 
onances. In addition, these polynomials will grow with 
energy and severely violate the unitarity bounds. 

It is, however, well known that the two caveats above 
can be fixed by the unitarization of the ChPT amplitudes 
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[14, 22, 25-28, 41-43]. In practice, one first unitarizes 
the partial waves of definite angular momentum J for 
the scattering of two pseudoscalars, which generates the 
poles associated to the a (980) for isospin I = 1 and 
/o(980) for isospin 1 = (the latter together with the 
/o(600), which is of little relevance in the context of this 
article). These unitarized scattering amplitudes are then 
matched to the form factors in such a way that the latter 
have the same poles as the scattering amplitudes, and 
therefore the same resonances. They do not grow with 
energy either and satisfy Watson's theorem of final-state 
interactions. 

Let us briefly review the general formalism before in- 
troducing the necessary modifications for the subsequent 
inclusion of isospin breaking. Assuming that only two- 
body intermediate states are relevant for the process, 
the unitarity conditions for the T-matrix, once projected 
onto partial waves of definite angular momentum J (see 
Appendix D), read 



ImT(s) = T(s)Z(s)T*(s) , 
Imr(i) = T(s)E(a)r*(s) , 



(7) 
(8) 



where 



T(s) = 



S(s) = — 
V ' 16tt 



Tn (a) T 12 (s) Tis(s) 
T 12 (s) T 22 (s) T 23 (s) 
Ti 3 (s) T 23 (s) T 33 (s) 

o-i (s) 
a 2 (s) 
o 3 (s) 



T(s) = 




(9) 



and Tij(s) is the partial wave T-matrix with angular mo- 
mentum J of the scattering between states i and j. The 
(Ji(s) = 2ki/yfs account for the phase space of the inter- 
mediate two-body states, where fcj is the center-of-mass 
momentum of each physically accessible state i. 

For illustration we have explicitly written the equa- 
tions for three coupled states i, j = 1,2,3, because in the 
problem at hand, the 1 = and J = isospin conserving 
process corresponds to this formalism with the identifi- 
cation 1 = 7T7T, 2 = KK, 3 = 7777, and the 1^ are just the 
form factors defined in Eq. (6) above. Nevertheless, in 
Refs. [23, 24] the calculations were performed in a two- 
channel approach neglecting the final state effects of -qrj 
rescattering, which will be included here for completeness 
(and consistency). Let us remark, finally, that the above 
formulae as well as the following results in this section 
are ready for a straightforward generalization to the case 
when the isospin states mix. In particular we will have to 
deal with six coupled channels once we introduce isospin 
violation in the next sections. 

As explained above, the ChPT partial wave T-matrices 
and form factors cannot satisfy Eqs. (7) and (8). Let us 
remark that Eq. (7) implies that T(s)" 1 = ReT _1 (s) - 
iE(s), so that, indeed, we only have to find an approxi- 
mation to ReT -1 . Actually, partial waves satisfying the 
coupled channel unitarity constraint are obtained from 



the following expression [14, 22]: 

T(s) = [I- T (2) (s) G(s)] _1 T (2) (s) , (10) 

where I is the identity matrix, is the leading-order 
ChPT T-matrix, T = T< 2 ) + 0(p 4 ), and G(s) is a diag- 
onal matrix whose elements Gi(s) are one-loop integrals 
corresponding to the two mesons of the state i propa- 
gating in the loop; detailed expressions are provided in 
Appendix B. Note that, in the physical (or charge) basis, 
G(s) is a matrix whose diagonal elements Gi(s) are an- 
alytic functions except for a right cut starting at each i 
state threshold. The imaginary part of Gi(s) is precisely 
oi (s)/16ir. Moreover, if T is reexpanded, one recovers the 
leading-order ChPT result including the correct imagi- 
nary part obtained at one loop. In addition, the factor 
[I — T^ G(s)]^ 1 generates the required poles associated 
with resonances. Alternative derivations of this unita- 
rization formalism make use of Lippmann-Schwingcr-likc 
equations [14] or dispersive approaches [25-27]. 

Unitarization thus provides the summation of the two- 
meson s-channel loops, since the G(s) functions yield the 
correct imaginary part and a cutoff that can be fixed 
to approximate the real part of ReT -1 , effectively ab- 
sorbing higher-order terms. This has been shown to 
be sufficient to reproduce the available scattering data 
on the scalar channel and generates the observed res- 
onances [14, 22]. Of course, with one natural cutoff 
this approximation is not always valid over the whole 
light resonance regime, and indeed next-to-leading-order 
terms are necessary to generate other states, like vectors, 
from unitarization [25-27]. 

Equation (10) shows how to unitarize scattering am- 
plitudes, but starting from this expression it is straight- 
forward to also obtain scalar form factors that satisfy 
Eq. (8), writing [23] 



T( S )=[l-T^( S )G(s)]- 1 R(s) 



(11) 



where R(s) is a vector of real functions free from any 
singularity, which can be determined from a matching to 
the next-to-leading-order ChPT calculation of the form 
factors T = rW +r( 2 ) +0(p 4 ) [44]. By reexpanding [I - 
T<- 2 \s)G{s)}- 1 ~ I + TW(s)G(s) + one can extract 
R(s) from 

T(a) = [l + T^G( S )]R(s) + 0(p i ) . (12) 

Since the G(s) integrals do have a residual cutoff q ma x 
dependence, not present in the renormalized next-to- 
leading-order ChPT calculation of T, this matching has 
to be performed at a given renormalization scale fi = 
l.2q max (see [22] for the relation between the cutoff and 
dimensional regularization of the G(s) functions). We 
provide the explicit expressions for R(s) in Appendix C. 
We will now discuss how this formalism needs to be ex- 
tended to include the effects of isospin violation. 
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III. ISOSPIN VIOLATION, THEORETICAL 
FRAMEWORK 

Let us focus on the production reaction J/ty — ► <f>ir°ri, 
and calculate the matrix element 

M%> = C (O|ss + A^m|A) . (13) 

In analogy to what we described above for the isospin- 
conscrving case, the matrix element will be written in 
terms of scalar form factors 

-V2BT^(s) = (0|fm|A> , 

(14) 

-V2BT s mi ( S ) = (ON A) , 

which obviously vanish in the isospin limit, as the source 
terms are isoscalar, while the final state has 7=1. We 
will start with a variant of the recent analysis of the reac- 
tion J/'fy — > (/>7r + 7r~ [24]. There the data reported by the 
BES collaboration [39] was analyzed using the same uni- 
tarized chiral approach, sketched in the previous section, 
for the meson-meson final-state interaction. Thus all that 
needs to be done to investigate the effect of isospin vio- 
lation is to replace the matrix G(s), the vector R(s), and 
the meson-meson scattering matrix used there by those 
including isospin violation. 

A. Form factor unitarization with isospin violation 

Unitarization does not actually rely on isospin conser- 
vation, but is just a formalism derived for partial waves 



that could have been applied equally well in the physi- 
cal (charge) instead of the isospin basis, which is indeed 
the most appropriate, once we allow for isospin breaking. 
Isospin-breaking effects in unitarized chiral effective the- 
ories have been studied before in the context of r\ and 7/ 
decays, see Refs. [45, 46]. 

The enhancement of isospin violation we are interested 
in is due to the fact that we are looking at the region 
of the two KK thresholds, where the dynamics is dom- 
inated by the resonances already generated within the 
chiral unitary approach. Isospin breaking is a small cor- 
rection to the isospin-conserving formalism, which simply 
amounts to increasing the number of distinct states and 
to slightly modifying the structure of the vertices. 

In the following sections we will thus perform the cal- 
culation of the different pieces needed in order to include 
isospin violation in the chiral unitary approach to the 
scalar form factors, as described in Sec. II B. In partic- 
ular, in addition to the three 1 = states, we now also 
have to consider two 1=1 and one 1 = 2 states. Hence, 
the T( 2 ) matrix is now six-dimensional, and its elements 
are easily obtained from the partial waves, shown in Ap- 
pendices A and D. The matrix of loop functions is also 
six-dimensional, as seen in Eq. (22). 

We can obtain all unitarized form factors for J/\& de- 
cays into a (j) plus two S'-wave pseudoscalar mesons in the 
isospin basis from the following equation: 



r£°(«) 
r£=°(s) 



( 



[l-TW(s)G(s)} 



\ 



R s w (s) + \<t>K(s) 
R°(s) + V$(«) 
R* K (s) + X^s) 
MR u K d {s) 
\^W + W«) + A*<(»)/ 



(15) 



where all the R functions can be found in Appendix C 
and in Eqs. (25) and (36). The 1=1 polynomial form 
factor terms oc A^, as well as the appearance of I = 
form factors ii™' s (s) in the nrj channel will be discussed 
in detail in Sec. IIIF. We neglect an isospin- violating 
production vertex with 1 = 2 since it would be sublead- 
ing in rrid — m u . Similarly, contributions oc \<i>R^ K (s) 
could be added to the 1 = production terms, but are 
neglected as second order in isospin violation. Let us 
remark that if we turn off isospin violation in and 
G(s) and set A^ = 0, R™ v (s) = R^ n (s) = 0, we recover 
the three coupled channel isospin-conserving case. 



For the study of a — f mixing, we are interested in 
the 1=1 form factors and more precisely in T^ 1 (s) 
in the region around KK threshold. Note that Eq. (15) 
is very convenient in order to switch on and off the dif- 
ferent contributions to isospin violation, and study their 
sizes, since all the isospin violation in vertices appears 
in T( 2 \ and the difference between charged and neutral 
loops appears in G(s). This will be studied in Sec. IV. 
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Why a — f mixing should dominate isospin 
violation 



In Ref. [20, 21] it was claimed that, if the invariant 
mass of the outgoing two-meson pair is close to the nom- 
inal mass of both / and do, then the mixing of the two 
should be, by far, the dominant isospin-violating effect. 
The argument was based on the fact that the two scalar 
resonances of interest are narrow and overlap and there- 
fore the effect of isospin violation as it occurs in the prop- 
agation of the scalar mesons is enhanced compared to 
mixing in the production operator. In Rcfs. [20, 21] the 
reasoning was presented for NN induced production of 
the scalar mesons. Here we adopt it for J/^S decays. For 
the sake of this argument only, we introduce the notion 
of explicit resonance (ao, /o) propagators; we emphasize, 
though, that such objects never appear in this form in 
our unitarized chiral amplitudes, where the correspond- 
ing poles are generated dynamically. 

We expect the effect of isospin violation in the produc- 
tion operator (oc in Eq. (15)) to scale at most as (to^ — 
m u )/m a . This is then followed by an isospin-conserving 
final-state interaction proportional to W ao ^TT V G ao , where 



G ao denotes the ao propagator and W a 



the ao de- 



cay matrix element — if we assume a scalar coupling of 
the ao to irrj, the vertex function W ao ^nri reduces to the 
effective coupling constant g aa iT-q (c.f. the corresponding 
couplings to the kaon channels defined in Eq. (1)). On 
the other hand, ao — /o mixing occurs in the propaga- 
tion, here parameterized by the various isospin-violating 
scalar form factors T™^. We may therefore use 



Wag-^ir,) G ao A Gf 



(16) 



Here A denotes the mixing matrix element, see Eq. (4). 
As it has been argue d in the introduction, A scales as 
^K\f{ m d — m u )/m s . For the / propagator very close 
to the KK threshold, we may use 



G fo 



1 



I, f 
./<) 



im fo r fo 



(17) 



which reduces to — i/(rrif a Tf a ) for s ~ m aa ~ TO / - 
the ratio of the two effects, and therefore an estimate of 
the theoreti cal uncertainty of the investigation, is given 
by m/ o r f a y/(md — m u )/m s /M] i , which is of the order of 
4% for the amplitude (for this estimate we used 50 MeV 
for the width). This estimate can now be tested within 
a dynamical approach. 



C. Lagrangians, mass splittings etc. 

We use the leading order chiral Lagrangian £( 2 ) includ- 
ing isospin-breaking/electromagnetic effects [47]: 



£(2) = —(D^UD^U^+xU^+Ux^+CiQUQU^) . (18) 



U collects the Goldstonc boson fields in the usual 
manner, F is the common pseudoscalar meson de- 
cay constant. The covariant derivative in particu- 
lar contains the coupling to photons, D^U = d^U — 
iQlAftjU] + and Q is the quark charge matrix, 
Q = ediag(2, — 1, — 1)/3. The field \ collects the quark 
masses, \ — 2B diag(m„, m,d, m s ) + We choose to ex- 
press isospin-breaking effects due to the light quark mass 
difference in terms of the leading order n r] mixing angle 



V3 m d - m u , 3 s 
e = — — +0((m d -m u Y) 



m s — m 



(19) 



Electromagnetic contributions oc C in Eq. (18) can be 
reexpressed in terms of the charged-to-neutral pion mass 
difference, 



(20) 



where the tiny strong mass difference cx (m<2 — m u ) 2 is 
neglected. Because of Dashen's theorem, the charged- 
to-neutral kaon mass difference can then be written at 
leading order as 



A 



K 



M 



K+ 



M 



4e 



K° 



A n -^={M 2 K -M*) . (21) 



As outlined in the introduction, the meson mass differ- 
ences, especially those of the kaons, naturally introduce 
a striking isospin-violating effect. Note that, as a con- 
sequence of the previous arguments, the matrix of loop 
functions G(s) is now six-dimensional, i.e. is diagonal in 
the charge basis 

G = diag(G 7r + 7r - , G^-o^o, G mi G K + K - , G K o^o, G^o^) , 

but has only a block diagonal form in isospin basis (see 
Appendix A). 

We have pointed out in Eq. (1) that the mass difference 
in the kaon propagators generates the leading isospin- 
breaking contribution in the ao//o resonance region. As 
we are about to calculate subleading isospin-violating ef- 
fects, one may wonder how accurately the unitarity cut 
contribution is described by the (leading-order) chiral 
unitary approach, and how much this description would 
change if we considered unitarized chiral p 4 amplitudes 
(see Refs. [25-27]). In this formalism we do not introduce 
the couplings of the ao and /o to different isospin chan- 
nels explicitly, which, as seen in Eq. (1), determine the 
strength of the leading part of the mixing amplitude. Ac- 
tually, within this approach, those couplings correspond 
to the residues of the poles that have been generated 
dynamically when fitting the data. In this sense the cou- 
plings are determined by the set of data that has been fit- 
ted to obtain the isospin-conserving part. The formalism 
as presented is very general and for the S'-waves, the qual- 
ity of such fits depends mainly on the data considered, 
and hardly changes with the order of the unitarized chiral 
amplitudes. Of course, whenever new data appears for 



the isospin-conserving reactions, the corresponding parts 
can be refitted and therefore allow for an improvement 
in the accuracy also for the isospin- violating amplitudes. 
However, the relevant observation is that at any time the 
possible effect of isospin violation in the couplings, as 
discussed in the following section, should be considered. 





FIG. 4: Photon exchange diagram that generates a singu- 
lar behavior at the K + K~ threshold (right). The full lines 
denote charged kaons, dashed lines arbitrary other mesons, 
and the wiggly line the exchanged photon. This is added to 
the standard K + K~ one-loop function (left) in the iterated 
bubble sum. 



D. Isospin violation in vertices 



E. Coulomb corrections 



Isospin violation in the Lagrangian in Eq. (18) does 
not only induce charged-to-neutral pion and kaon mass 
differences, but also affects the (tree-level) scattering 
amplitudes, which receive corrections to their isospin- 
conserving expressions. This formalism allows one to 
calculate isospin violation in scattering consistently with 
the analysis of the masses. A well-known example for 
the importance of such effects is the sizeable correction 
in relating 7T7t scattering amplitudes at threshold to the 
isospin scattering lengths [48]; similar corrections have 
also been calculated e.g. in irK scattering [49-52]. In 
addition, some transitions only take place at all in the 
presence of isospin violation, the most prominent being 
the decay -q ^ 3n [53, 54]. 

The isospin-violating scattering T-matrix has to be cal- 
culated in the particle or charge basis; the complete list 
of amplitudes, linking all charge- and strangeness-neutral 
channels, calculated at leading order within ChPT, i.e., 
is given in Appendix D. Note that we 
now have six coupled states, but as we have commented 
in Sec. II B, the same unitarization formalism applies. 
Despite obtaining our calculations in the charge basis, it 
is still convenient to recast them in the isospin basis, as 
we assume (for the moment) production of a pure 1 = 
state, while 7r°?7 is I — 1. We show its relation to the 
charge basis in Appendix A. Obviously, the matrix of 
scattering amplitudes in the isospin basis is not block di- 
agonal with respect to I = 0, 1, 2 anymore, but allows 
for transitions between different isospin quantum num- 
bers; these transition matrix elements scale with either e 

9 

or e . 

In the resonance picture of Eqs. (1) and (3), isospin vi- 
olation in the meson-meson vertices induces the isospin- 
violating resonance couplings 9 RK ^- As we have no 
strict counting scheme for energies in the resonance re- 
gion, these resonance couplings are only modeled this 
way, and may receive corrections from higher orders, al- 
though one should take into account that unitarization 
is necessarily taming their effect. They are not fixed by 
isospin-symmetric data in the way the cut contribution 
due to kaon mass differences are, and therefore have to be 
pinned down directly from isospin-violating decays, like 
the one discussed here. 



So far, all interaction terms derived from Eq. (18), en- 
tering the matrix in the unitarized final-state in- 
teraction, are pointlike four-meson vertices. We note, 
however, that Eq. (18) generates another tree-level di- 
agram contributing to meson-meson scattering, namely 
one-photon exchange between charged mesons. It is ob- 
vious that this nonlocal interaction cannot be taken into 
account on quite the same footing. 

As all initial- and final-state particles in the decay 
J/ty — ► <j)ir°rj are electrically neutral, photons can only 
enter inside charged-meson loops. The diagram shown in 
Fig. 4 is the only one at 0(a) that is enhanced at the 
K + K~ threshold; we neglect all other, nonenhanced di- 
agrams. Our prescription is then to replace the charged- 
kaon loop function Gk+k- (s) m the unitarization sum 
by the sum of this and the one-photon exchange graph, 



G 



K+K 



{s)~,G K+K -{ 3 ) + G 1 K ^{ s ) 



(23) 



In the threshold region, the exact expression for the one- 
photon-exchange diagram [55] can well be approximated 
by the threshold-expanded form [56] , which reads 



G 



-4M|. + | 



log 2 



21C(3) 
2tt 2 



m9( S - 4M| + )| + 0((s - 4A4 + ) 1/2 ) 



(24) 

A justification for this being the leading threshold behav- 
ior can also be given in the framework of a nonrelativistic 
theory [57]. 

We neglect exchange of multiple photons inside the me- 
son bubble, although these are in principle also enhanced 
close to threshold. A resummation of multiphoton ex- 
change is necessary once the parameter 



4M| + 



-1/2 



is not small anymore; it, e.g., becomes as large as 0.1 
for y/s = 2M K + ± 0.7 MeV. The size of the more-than- 
one-photon exchange amplitude relative to one-photon 
one stays below 5% outside a window of ±5 MeV around 
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the charged kaon threshold. As the energy resolution for 
potential experiments is expected to be of this order at 
best, neglecting higher-order photon graphs seems well 
justified. This is in line with the findings of Ref. [58]. 

Note that we only take photon exchange into account 
for the charged kaon loop graphs. We have checked 
that the corresponding modification inside the charged 
pion loops leads to no visible modification in the phys- 
ical region of the process under investigation, i.e. above 
ir°r] threshold, and in particular not in the energy region 
around KK threshold considered here. 



J/* 

\ 



X 



~^{uu — dd) 

FIG. 5: Vector-meson mixing contribution to A^. The box 
denotes an (isospin-conserving) coupling of J/^ to p° and 
the 1 = 1 scalar source term, the cross the isospin-violating 
p° — (j> mixing term. 



F. Isospin-violating production operators 

It has been argued in Sec. Ill B that isospin violation 
in the final-state interaction ought to be the dominant 
effect for the production of the Tr°r] final state in the en- 
ergy region around the KK threshold. We can check this 
assumption explicitly by allowing for isospin breaking in 
the production operator. This occurs in two forms: due 
to mixing, the ir°r] final state has an / = component, 
i.e. nonvanishing form factors with the nn and ss cur- 
rents that scale with the mixing angle e exist already at 
tree level; and we may allow for an additional 1=1 
component in the scalar source terms given in Eq. (5). 



1. Isospin-violating scalar form factors at tree level 

The mixing of ir° and r\ leads to nonvanishing form 
factors r™^(s) already at tree level. At leading order in 
the chiral expansion, where the propagators for ir° and 
r\ can be diagonalized with a single mixing angle e, the 
matrix element (O\uu+dd+ss\ir r]} has to vanish, leading 
to the relation y/2T™ v (s) + T s „ v (s) = 0. We find 

Kv(*) = -^ R ^( s ) = ^ + 0(e 3 ), (25) 

as r(s) = R(s) + 0(p 2 ), see Eq. (12). These components 
are added to the vector of production operators R of 
Eq. (15). In contrast to the 1 = 1 scalar source term dis- 
cussed below, this effect does not induce any additional 
uncertainty, but comes with a fixed coefficient. 

In principle, the task of this investigation consists pre- 
cisely in the determination of the form factors r™ (s), 
(s), and T"^(s). In contrast to what was done in 
the isospin-conserving case, however, we only match the 
ChPT expressions for these form factors at leading order. 



2. 1=1 scalar source term 

We first introduce an I = 1 scalar source, generalizing 
the Lagrangian in Eq. (5) to 



C = 



ss + X^nn H — j= (uu — dd) 
v2 



(26) 



In order to obtain an order-of-magnitude estimate of 
the size of the new parameter A<^, we invoke a version 
of vector-meson dominance: we assume the dominant 
(isospin-conserving) production of a p° in association 
with the operator (uu — dd)/y/2, with subsequent p° — </> 
mixing; see Fig. 5. Following Ref. [24], we write the in- 
teraction Lagrangian in the form 



g^{(V,S) + ^(V,)(S)} , (27) 



where collects the vector meson fields, S is the matrix 
of scalar sources, and S = S — (S)/3. Considering just 
the flavor-neutral vector mesons, we can rewrite Eq. (27) 
as 



L = 



* M |C0 4>f,(ss + X^nn) + w M (ss + X^nn) 



(28) 



where C u , X u , C p can be expressed in terms of Cp, A</> 
according to 

Cuj = X<f,C<f, , A w = ' C P = (l — —^jCp . (29) 

p — <f> mixing is actually assumed to proceed via sub- 
sequent p° — to and lo — 4> mixing. Neglecting the finite 
widths of the vector mesons, we find a coupling of the <f> 
to the (uu — dd) / \J~2 operator in terms of the p° — lo and 
u — 4> mixing angles 6^, 0^ of the form 



Cf^Xft 

and therefore 

Xih = 



(Ml-Ml)(Ml-Ml) 



(Ml-Ml)(Ml-Ml)\ V2J 



(30) 



(31) 



Plugging in the central values for the mixing angles from 



Ref. 



e 



f)uJ 



-3.75 x 1CT 3 GeV 2 , 9, 



1(T 3 GeV 2 , we obtain 

L, w - 0.5 x 10~ 3 



25.34 x 



(32) 
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where the uncertainty due to errors in the input numbers 
is about 20%. It is interesting to note that a correspond- 
ing estimate based onw-^ mixing of the A^ parameter, 
which violates the Okubo-Zweig-Iizuka rule, would lead 
to 



A^ = 



e 



U)(j> 



M 2 - M 2 



0.06 



(33) 



which is only about a factor of 2 smaller than the fit 
result [24]. We therefore assume that the above order- 
of-magnitude estimate of A^ should be comparably ac- 
curate, and consider it a very conservative estimate to 
vary the strength of the isospin- violating production op- 
erator within a range for A^ increased by a factor of ±10 
compared to Eq. (32). 

In order to follow the formalism used earlier, we gen- 
eralize the matrix element Eq. (13) to 

MJ 1 = C^Olss + X^nn+^iuu-dd^rj) , (34) 

and define the additional I = 1 scalar form factors T^f, 
T"^ according to 



-V2BT u K d (s) = (0\^=(uu-dd)\KK) I=1 , 
v2 

-V2BTf v (s) = (0|^(««- dd)|A) J= i • 



(35) 



As we are only interested in an order-of-magnitude es- 
timate for the effects of the 1 = 1 production operator, 
we refrain from doing the complete one-loop calculation 
of these form factors and only unitarize the lowest-order 
(0(p 2 )) results for these, which we find to be 

Rg(8) = -j=, Rf ri (s) = -^= + 0(e 2 ). (36) 

We note that the Feynman-Hellman theorem implies the 
relation 



T) n ( s ) - Lrs ( s ) - ^ B ud ( s ) 



(37) 



at tree level, which can easily be checked to be fulfilled 
by Eqs. (25) and (36). 



G. Further possible background terms 



found to be much smaller than the signal and therefore 
will not be considered any further. 

The latter, however, turned out to be of the order of 
or even larger than the signal and will now be discussed 
briefly. In the unitarized chiral approach, only the light- 
est pseudoscalar mesons appear as dynamical degrees of 
freedom. The effects of all other mesons, e.g. the K* , are 
integrated out and are considered through local counter- 
terms. The corresponding expressions can be deduced 
from those with dynamical heavy particles by formally 
taking the infinite mass limit. 

According to Ref. [38] the transition in the decay chain 
J/\E' — > (K*K + h.c.) — > 4>clq appears through a trian- 
gle loop that contains two kaon propagators and one K* 
propagator; isospin violation then emerges through the 
kaon mass differences. In the effective theory description 
the vector meson propagator needs to be replaced by a 
point interaction — as a result the whole transition is to 
be regarded as part of the isospin-violating scalar form 
factor and therefore, from our point of view, as part of 
the signal. In this sense it appears natural that the cor- 
responding transition rate is of the order of the estimated 
signal. 



IV. RESULTS AND DISCUSSIONS 

The differential decay rate for the process J/* — > cfy-K^rj 
is related to the absolute square of the matrix element 



Ml r] = -s/2BC 4> 



(38) 

with the form factors (s) T^(s), and T^(s) given by 
the matrix relation Eq. (15), by 



pol 



\M 



^, 8 ,m|)a( 8 ,^,^) 

d^7s _ 16^Mj /vp (27r) 3 

(39) 

where A(x, y, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz) is the 
usual Kallen function, and i^pd is the kinematical factor 
that takes the average and sum over polarization states 
of J/* and into account, 



p,p' 



i + 



(M 2 ^ + M 2 - S y 

8Mj /iS M 2 



(40) 



In Ref. [38], two candidates for background terms to 
an a — f mixing description of the decay J/ty ~ * (p-K^-q 
were estimated, namely J/^f —* 7* — > cf)7r r) and J/^f —* 
(K*K + h.c.) -» W 

The former is clearly not covered 
by our description of this decay in terms of scalar form 
factors and would necessitate a generalization of the pro- 
duction mechanism beyond qq operators; luckily it was 



The data for J/* -> ^tt+tt" [39] analyzed in Ref. [24] 
do not provide normalized differential decay widths, but 
only event distributions, dN/d^fs cx dT/d^/s, with an 
unknown constant of proportionality. Accordingly, we 
cannot predict a normalized differential decay width for 
J/ty — ► (f>ir°r) either, but only a relative width, normal- 
ized by the isospin-conserving 2-pion decay channel. We 
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choose to perform the normalization according to 

V norm 



(41) 



2M 2 K+ =F 25 MeV cov- 



where the energy range W\.2 
ers the peak region of the clq — fo mixing signal. The 
numerical input on masses and coupling constants that 
enter our calculations are given in Appendix E. We re- 
mark that the coupling constant as well as the form 
factor normalization constant B, see Eq. (38), cancel out 
in the ratio in Eq. (41) and do not have to be specified 
in our analysis. We wish to emphasize that there are 
no new fit parameters in this analysis: the unitarization 
procedure in the chiral unitary approach at this order 
contains one single parameter, the cutoff g max , which has 
been adjusted as to reproduce the S'-wave meson-meson 
scattering data in all physical channels (i.e., in particular 
the resonance pole positions) as well as possible. All fur- 
ther parameters, both the parameter in Eq. (5) as well 
as the 0(p 4 ) chiral low-energy constants L\ in the vari- 
ous scalar form factors, only serve as a parameterization 
to describe the isospin-conserving J/\I> — > <f)irn data. 

Of course, in order to predict realistic relative count 
rates, one would have to take into account different de- 
tector efficiencies for the two different final states, which 
we do not do here; it will be straightforward to imple- 
ment them, once such experimental specifications become 
available. 

As stated above, the tt T] invariant mass distribution of 
the reaction J/^f — > (j}TT°r] is proportional to the absolute 
value squared of the mixing matrix element, at least as 
long as we do not include any isospin violation in the 
production operator. In Fig. 6 we show our prediction 
for the TT°r] invariant mass distribution when different 
isospin-violating effects are included in the propagation 
of the scalar mesons. To produce the dashed line only the 
kaon mass difference was included. The resulting curve 
is in qualitative agreement with those of Refs. [18, 19]. 
Including isospin violation in the vertices, we find the 
dot-dashed curve; adding one-photon exchange according 
to Sec. HIE leads to the full result, given by the solid 
line. We find that the effect of photon-exchange is rather 
small even in the threshold region, and certainly smaller 
than the modifications due to isospin violation in the 
meson-meson scattering vertices (for a more detail study 
of photon effects see Ref. [58]). The signal for a — f 
mixing is therefore significantly enhanced compared to 
the original estimate given in Ref. [18]. 

Next we investigate the possible effect of isospin vi- 
olation in the production operator. As one can see in 
Fig. 7, in the region around the kaon threshold the in- 
variant mass distribution in the -K°r] channel is by far 
dominated by the isospin violation in the propagation. 
The admixture of an isospin-violating production oper- 
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FIG. 6: Predictions for the normalized J/\I/ — » differen- 
tial count rate per 10 MeV, relative to 1000 J/* — > 4>it + -k~ 
events in the peak region 2M^. + ± 25 MeV as denned in 
Eq. (41). We show the successive inclusion of the different 
isospin-violating effects in the final-state interaction. The 
kaon mass difference alone leads to the dashed line, additional 
inclusion of isospin violation in the strong vertices produces 
the dot-dashed curve. As the full result, adding one-photon- 
exchanges in the bubble sum, we obtain the solid curve. The 
top panel shows the most relevant region near the two KK 
thresholds; in order to make the enhancement in this region 
more obvious, the bottom panel shows the whole energy range 
from threshold up to 1.2 GeV on a logarithmic scale. 



ator as estimated in Sec. IIIF actually produces only a 
rather narrow band very close to the result without such 
an operator. This nicely confirms the corresponding es- 
timates provided in Ref. [20, 21]. 

In order to demonstrate the enhancement of isospin- 
breaking in the two-kaon threshold region, as well as due 
to different mechanisms, we show the predicted event 
numbers in Table I, normalized to 1000 J '/\& — > (/>7r + 7r~ 
events in the peak region, for three distinct kinematical 
regions: for the two-kaon threshold region of total width 
50 MeV, below down to threshold, and above up to 
1.2 GeV. We want to point out that the signal around 
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FIG. 7: Predictions for the J/ty — > <f>-K°ri differential count 
rate per 10 MeV (normalization as in Fig. 6) , with (band) and 
without (full line) inclusion of an isospin-violating production 
operator (IVPO). The central curve is shifted by the isospin 
violation in the scalar form factors as discussed in Sec. IIIF 1, 
the band is due to the I — 1 scalar source with strength 
limited by A = ±5 x 1CT 3 , see Sec. IIIF 2. As in Fig. 6, the 
top panel shows the two-kaon threshold region, the bottom 
one the whole energy range from threshold up to 1.2 GeV on 
a logarithmic scale. 



1 GeV is enhanced by nearly 50% compared to the orig- 
inal estimate based on the kaon mass difference effect 
alone. The reduction due to photon exchange visible 
in Fig. 6 is about 10%, the corresponding event num- 
ber in Region II of Table I without photon exchange 
would be 33.6. Obviously, all mechanisms other than the 
kaon mass difference also lead to a large relative increase 
of mixing outside this central region, in particular the 
isospin-violating production operator; this is also clearly 
seen in the logarithmic plots over a wider energy range 
in Figs. 6 and 7. However, a relative enhancement of 
the signal by more than two orders of magnitude close 
to the two-kaon thresholds remains. Finally, we wish to 
point out that 29.0 J/^ — > 4>TT°r] events in the resonance 



TABLE I: Event estimates for J/$ — > <f>n°ri for different en- 
ergy regions and different isospin-breaking effects. The energy 
regions are: Region I: M w o + M v < ^/s < 2M K + - 25 MeV; 
Region II: 2M K+ -25 MeV < yfs < 2M K+ +25 MeV; Region 
III: 2M K+ + 25 MeV < y/s < 1.2 GeV. "AM K " refers to the 
original estimate [18], assuming the kaon mass difference as 
the only source of ao — fo mixing. "Full FSI" labels the model 
with all isospin-breaking effects in the final-state interaction 
(FSI) included, while "Including IVPO" also incorporates the 
isospin-violating production operators discussed in Sec. IIIF. 
All event numbers are relative to 1000 J/* — > 4>-k + -k~ events 
in Region II. 
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region relative to 1000 J/^ — > </>7r + 7r~ in the same kine- 
matic range may look like a mere 3% effect of isospin 
violation, but, due to the lack of interference, it corre- 
sponds to about 17% isospin breaking on the amplitude 
level. Therefore, oo — fo mixing leads indeed to a very 
sizeable isospin-violating signal. 

As a side remark, we briefly comment on the form fac- 
tors r™ J) (s), r* ?? (s) individually, which in principle may 
occur in a different relative combination when ao — fo 
mixing is investigated in the context of a different decay 
or production mechanism. In Figs. 8 and 9, we show the 
absolute values of the form factors as well as the differ- 
ence of the phase motions of the two (both including all 
isospin-breaking mechanisms discussed earlier). We find 
that the shape of (the absolute values of) both form fac- 
tors is very similar, and that the phase difference, while 
showing remnants of the cusps, varies rather mildly over 
the two-kaon threshold region (by altogether less than 
15°). We conclude from these observations that, if the 
relative strength parameter is replaced by a differ- 
ent, possibly complex (but still approximately constant) 
value, we would not expect the mixing signal to be wildly 
different from what we predict here. 



V. SUMMARY 

In this work we have improved the theoretical under- 
standing of the phenomenon of a — f mixing. We have 
confirmed that the dominant mixing effect comes from 
the kaon mass difference in line with Ref. [18]. It is 
therefore possible to extract independent information on 
the effective couplings of ao(980) and /o(980) to kaons, 
which contain important structure information [15], from 
the mixing matrix element. 

We have applied our formalism to the reaction J/ty — > 
4>ir°r) and give what we consider the best and most com- 
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FIG. 8: Absolute values of the scalar form factors r™ (s) 
(full line), r^(s) (dashed line). Note that, to facilitate the 
comparison, the latter has been scaled down by a factor of 2. 
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FIG. 9: Relative phase of non-strange and strange nr/ scalar 
form factors, Arg[r™„(s)] - Arg [T%{s)] . 

prehensive prediction to date for this decay channel. The 
corresponding measurement will soon be possible with 
the upgraded BESIII detector [38]. 

In addition we addressed the following items: 

1. As it is illustrated in Fig. 6, although the kaon 
loop effect dominates the isospin-violating signal, 
the presence of isospin-violating coupling constants 
introduces a significant additional enhancement by 
roughly 50% (corresponding to a 20% effect on the 
amplitude level). To the order we are working, 
isospin violation in the vertices and in the kaon 
masses are of the same origin (i.e. are calculated 
from the same chiral Lagrangian), and no addi- 
tional parameters enter. 



2. The effect of soft photon exchange in the meson 
propagation is small, even in the signal region, and 
amounts to a reduction of the signal by 10% (or 5% 
on the amplitude level). 

3. Neither of the two additional effects studied in 
this work, associated with an isospin-violating pro- 
duction operator, distorts the shape of the signal 
severely, c.f. gray band vs. solid line in Fig. 7. 

4. We have confirmed that the isospin violation in the 
final-state interaction, which can be identified with 
a — /o mixing, is much more important than isospin 
violation in the production operator. We therefore 
confirm the corresponding claim of Refs. [20, 21] for 
the particular case of the reaction J/^f — > §ir Q r\. 

Our final results are presented for a particular reac- 
tion. However, the formalism to calculate the propa- 
gating a — /o system in the presence of isospin viola- 
tion is very general and therefore could also be used for 
the analysis of other experiments, once data is available. 
Especially items 1-3 of the list given above are indepen- 
dent of the reaction studied, and ought to be distinctly 
different from predictions of models viewing ao and f 
as qq or four-quark systems, which allow for direct mix- 
ing without two-kaon-thrcshold enhancement. We hence 
conclude that it is indeed possible to measure and inter- 
pret a — fo mixing. 
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APPENDIX A: CHARGE VS. ISOSPIN BASIS 

The two-meson states of definite isospin (in the S- 
wave) are related to the states in the charge basis by 
the following relations: 
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/ |7T7r) /=2 ^ 
|7T7r)j =0 

\w)i=o 
\KR)i=o 
\KK) I=1 
\ |A)/=i ) 



J=0 




\TT°n°) 

\vv) 

K+K~) 
\K°K°) 

I A) J j 



(Al) 



=o 



where \ir + ir ) denotes the symmetrized combination of 
\tt + tt~) and |7r _ 7r + ). The additional factors of 1/V2 for 
the 7T7T and the r\r\ states account for Bose symmetry. 
Using the isospin basis for our calculation is particularly 
convenient as we consider production of states with def- 
inite isospin (7 = mostly), and the final state n°r] has 
definite isospin 1=1, too. 



APPENDIX B: THE LOOP FUNCTION G 

For two mesons A, B with masses M A and Mb prop- 
agating in the loop, the elementary two-point function is 
given by 



Gab(s) 




. i + n A m a 
log TTI^- log M^ 



v (i/Cl A + s - A)(z/Q B + s + A) 
2s ° S (vfl A - s + A)(t/ft B - s - A) 



(Bl) 



log 



Q 2 

lrnax 



M A M B 



[i + n A ) (i + n B ) 



where A = M% - M%, v 2 = s 2 - 2s(M 2 + Mf) + A 2 , 



&A,B = y 1 + -^l ;S /9max ) and lrnax IS the Cutoff. The 

form in Eq. (Bl) is correct and unambiguous below the 
pseudothreshold, s < (M A — M B ) 2 , and has to be con- 
tinued analytically into other kinematical regions. As 
described in the text, the determination of Gi(s) has to 
be chosen so that on the physical, right hand, cut of the 
i-th state, ImGi(s) = <Ji(s) /16tv = ki/(8iryfs), where ki 
stands for the corresponding center-of-mass momentum. 
In the equal-mass limit M A = M B = M, Eq. (Bl) re- 
duces to G A=B (s) = G(s), 



G(s) = 



1 



16tt 2 



clog 



crSl + 1 
aVL- 1 



2 log 



M 



-(1 + 0) 



where a = ^1 -4M 2 /s, Q = ^l + M 2 /^ 



APPENDIX C: POLYNOMIAL REMAINDERS OF 
THE SCALAR FORM FACTORS 

The polynomial remainders of the scalar form factors 
as defined in Eq. (12), up to 0(p 4 ) in ChPT, are given 



as 



K(s) = yf {l + ftrW + - LI) 

8(2M| + 3M 2 ) 



F 2 



■(2L5-L3 



4s 
F 2 



(2LI + LI) 



v/lt 3 l[K) F 2 ^ 8 ~ 5 > 
+ 8(eM| + M;) g 4^ | 



+ 8 ^ 2 + ^ (2^-^) + ii(2^ + ^), 



K 



16M 2 . . 8(10M 2 - M 2 ) , 

v 



F 2 
n 



128( m|-m;) (3L7+l;)+ 4, (6l; + £8 



i^(«) = -H|l + 2/i Jf 



16M„ 2 ,___ _ 4(8M| + M 2 ), orl 

17 - 1 rj 



F 2 



^|^(3. 7+ ^) + |(3^ + 2^)} ) 

(CI) 



(E52) wnere are tne tadpole loop functions defined as 



M 2 



log . 



32TT 2 FiFj ~ n 



(C2) 



and this particular choice of decay constants ifyj is in 
accordance with Ref. [60] . The L[ are the standard low- 
energy constants defined in Ref. [61]. 
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APPENDIX D: SCATTERING AMPLITUDES 
WITH ISOSPIN VIOLATION 

In this appendix we give a complete list of S- wave pro- 
jected four-meson amplitudes derived from Eq. (18) for 
those channels that are charge- and strangeness-neutral, 
i.e. reactions linking the channels (1) 7r + 7r~, (2) tt°tt°, 
(3) f]t], (4) K+K~, (5) K°K°, (6) tt ^ including isospin 
violation; see also Ref. [45]. Our normalization for the 
S'-wave projection is given by 

T^=°( S ) = \ J ^dzT ab {s,t(s,z),u(s,z)) . (Dl) 

a, b refer to the numbering of the possible in- and out- 
channels 1-6 as given above; of course one has T a b = Tb a . 
In the following list, for reasons of brevity we suppress 
the J = superscripts: 



Tn — 2T 14 — T 44 — 



s + 4A„ 
2F 2 



'12 



'■13 



F 2 

1 M 2 

,, 22 — J 23 — J 66 — -^pj ' 



Tib = 
Tie = 



T45 — 2^55 



s 

af 2 



3F 2 

T24 = T25 = 



(4M 2 3s) , 



2e 



(1 ± 2V3e) T ~^M 2 K , 



AF 2 



T26 — 



T 



36 



V3F 2 



133 



AM 2 - M 2 



3F 2 



T 35 
T46 
T56 



(1 + V3e) , 



AF 2 \ + V3~J 

V3^ 



4 V / 3F 2 V V3 



4V3.F 1 



^-rt 1 '^!- 



All amplitudes are normalized by a factor of 1 / F 2 which, 
at this accuracy, can arbitrarily be identified with any 
meson decay constant. For numerical evaluations we use 
the convention that the overall 1/.F 2 factor is replaced 
by one 1 / \J~F\j, factor for every external meson (j) in the 
process concerned; for a discussion on how and why this 
choice yields the best description of data see Refs. [28, 
62]. 



APPENDIX E: NUMERICAL INPUT 



In our calculations we use the masses M n + = 
139.57 MeV, M T o = 134.98 MeV, M K + = 493.68 MeV, 
M K o = 497.67 MeV, M v = 547.8 MeV, M J/qi = 
3097 MeV, M = 1020 MeV. In addition we need the de- 
cay constants F n = 92.4 MeV, F K = 1.22F„ and F v = 
1.3F V . The leading order n°r] mixing angle is e = 0.01. 
For the low-energy constants of order p 4 [61] needed for 
the polynomial terms of the scalar form factors, we use 
the numerical values L\ = 0.84 x 10~ 3 , L\ = 0.52 x 10~ 3 , 
LI = -0.18 x 10~ 3 , L 7 = -0.40 x 10~ 3 , U 8 = 0.15 x 10~ 3 
(all given at a scale \i — M p ), as obtained in a three- 
channel generalization [60] of the formalism presented in 
Ref. [24]. From the same fit, we use the relative strength 
parameter of non-strange to strange scalar source terms, 
A0 = 0.117. In the loop function G(s) (see Appendix B), 
we employ a cutoff q ma x — 0.95 GeV. 
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